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We establish a set of nonequilibrium quantum phase transitions in the Dicke model by considering 
a monochromatic nonadiabatic modulation of the atom-field coupling. For weak driving the system 
exhibits a set of sidebands which allow the circumvention of the no-go theorem which otherwise 
forbids the occurence of superradiant phase transitions. At strong driving we show that the sys- 
tem exhibits a rich multistable structure and exhibits both first- and second-order nonequilibrium 
quantum phase transitions. 
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The Dicke model (DM) is a paradigm of collective be- 
havior in quantum mechanics [l| and describes the inter- 
action of N two-level atoms with a single-mode bosonic 
field. The DM undergoes a quantum phase transition 
(QPT) at a critical atom-field coupling, the same order 
of magnitude as the atomic level splitting Q • Whereas 
in cavity QED it is forbidden by a "no-go theorem" [3, 0] , 
the superradiant QPT has recently been realized experi- 
mentally in a simulation of the DM with a Bose-Einstein 
condensate in an optical cavity @ . 

Even under an adiabatic evolution many-body systems 
can exhibit highly nontrivial behavior [JQ- In particu- 
lar, spontaneous symmetry breaking at the DM QPT in 
a Bose-Einstein condensate was observed experimentally 
by adiabatically crossing the critical point [l3|. Fur- 
thermore, nonadiabatic evolution of quantum systems 
has attracted considerable interest, both theoretical and 
experimental [TTI-[T^ . 

In this Letter we study a driven version of the DM in 
which we assume a time-dependent atom-field coupling. 
In comparison with previous works which consider nona- 
diabatic modulation of a single two-level system Il5l4l8l 
or an TV-atom system under adiabatic modulation [7|, ll9[ , 
we address here the fundamental issue of the influence of 
nonadiabatic modulation on QPTs. 

In the limit of weak driving strength we show that 
our driven DM exhibits a set of new nonequilibrium 
normal-superradiant QPTs (QPT sidebands) when the 
driving is near resonance with the excitation energies 
of the undriven system. The QPTs are of second 
order and similar in kind to the original Hepp-Lieb 
transition [2[. We show that these nonequilibrium 
QPTs are not forbidden by the no-go theorem, thus 
bringing the otherwise-forbidden Dicke-type QPT in 
the realm of observability in cavity and circuit QED 
setups. Our analysis also allows us to go beyond this 
perturbative regime and investigate the limit of strong 
driving. In this regime we show that, in comparison 
with previous proposals for driven QPTs ^ 22^ , the 
nature of criticality changes dramatically. We find a 
rich nonequilibrium phase diagram replete with a host 
of macroscopically-distinct meta-stable phases and a 
nonequilibrium first-order QPT with no analogues in 
the static case. 



Driven Dicke model. — Following Dicke we de- 
scribe an ensemble of N identical, distinguishable two- 
level atoms (level splitting wq) by means of collective op- 
erators Ja]a € }, which obey the angular mo- 
mentum commutation relations (with cooperation num- 
ber j = N/2). These atoms interact with a bosonic 
mode of frequency ui via a dipole interaction. With the 
atom- field coupling strength time-dependent, we obtain 
the driven DM: 

Hit) = Loa^a + uJoJ. + ^ia'' +a){J+ + J-), (1) 

V-/V 

and in the following we shall consider a monochro- 
matic modulation with a static contribution: g(t) = 
g + Ag cos fit. 

Normal-phase stability. — In the thermodynamic limit, 
N oo, the undriven (Ag = 0) DM exhibits a QPT at 
a critical coupling gc = y/ZjZjo/2, where the ground state 
changes from an unexcited normal phase to a symmetry- 
broken superradiant phase in which both the field and 
atomic collection acquire macroscopic occupations 0, Q . 
We begin our analysis by investigating the stability of 
this normal phase under driving. To this end, we con- 
struct a normal-phase effective Hamiltonian in the same 
way as in Ref. [3| for the undriven case: we make a 
Holstein-PrimakofF representation of the angular momen- 
tum algebra in terms of bosonic operators b, b^ and take 
the thermodynamic limit, assuming b N'^ . The re- 
sult is the driven normal-phase Hamiltonian H]S!p{t) = 
LJob^b + uja^a + g{t){a'' +a){b'' +b), which describes fluc- 
tuations about the vacuum cavity state and unexcited 
atomic ensemble. In the Heisenberg picture, the equa- 
tions of motion for the normal coordinate operators of 
this model, q±{t), read [11] 



± 2ujAgcosQt 



q±{t) = Q, 



(2) 



with e± = y/iu^ zt 2guj, the excitation energies of the un- 
driven normal phase and where we have set = lo for 
simplicity. Equation ^ represents two uncoupled Math- 
ieu equations [2J|. In classical dynamics, the Mathieu 
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FIG. 1: Stability diagram of the nonequilibrium normal 
phase on resonance uj — ojq for (a) weak driving, — 0.15, 

and (b) strong driving, Ag/Q, — 0.4. The colored zones cor- 
respond to stability, white zones, instability. In the undriven 
case, the separatrix between stable and unstable zones is given 
by Qc — ui/2 on resonance (dotted red line). For weak driv- 
ing instability zones open up around the resonance conditions 
e± = k{^}/2) for fc = 1,2,3 (dotted black lines), which grow 
and begin to dominate for large driving. 



equation exhibits the phenomenon of parametric reso- 
nance and has stable and unstable solutions whose loca- 
tion is given by the Arnold tongues [ISl- The manifesta- 
tion of parametric resonance in the quantum regime has 
also been studied e.g., j25| . 

In our case, when both normal modes q± are stable, 
Hamiltonian _ffNp(i) permits bound solutions, localized 
around unoccupied field and atomic modes. Either scaled 
field occupation, 2{a^a)/N , or scaled atomic occupation, 
2{b''b)/N = 2{Jz)/N + 1, may be taken as the order pa- 
rameter for this system, and here they are both zero 
(A^ — >■ c»). When unstable, Hj<sp{t) possesses only un- 
bounded solutions and ceases to be a valid approxima- 
tion to the full Hamiltonian (admitting the possibility 
that the order-parameter becomes finite) . 

Figure [1] depicts these stable zones as a function of 
the static parameters g and uj — the colored zones cor- 
respond to stability, white zones, instability. Without 
driving, Ag = 0, the stability-instability transition corre- 
sponds to the standard DM QPT along the line g = oj/2 
(on resonance). Increasing the driving strength Ag has 
two effects: it leads to a shift of this critical line — as a 
consequence of parametric stabilization |24| — and more 



normal phase. The precise locations of these zones can be 
obtained from the known behavior of the Arnold tongues 
[23 |. When they first appear, these zones are located 
around resonance between driving frequency and (un- 
driven) excitation energies: kQ = 2s± with integer fc > 0. 
A similar instability for fc = 1 was briefly discussed in 
the dispersive limit u) ^ ujq in Ref. [26j . For sufficiently 
small driving, the width of the fcth unstable zone scales 
like (2/17)2''^-i(A5 uj)''. For large Ag the new instability 
zones dominate the parameter space [Fig. [TJb)]. Just as 
the change in stability of the undriven normal-phase ef- 
fective Hamiltonian indicates the occurrence of a QPT, 
we interpret the change in the stability of i?NP (t) as ush- 
ering the occurrence of a nonequilibrium QPT. 

Nonequilibrium QPT sidebands. — To obtain more in- 
formation about the nonequilibrium QPTs and the un- 
stable zones of Fig. [1] we employ Floquet theory [23] and 
make use of a generalized rotating wave approximation, 
similar to that found in Refs. [H, [l^]. Motivated 
by the fact that for small static coupling g, the fcth in- 
stability zone arises close to cj « kfl/2, we perform a 
canonical transformation of Eq. ([T]), H{t) — >■ Hk{t) — 
Ul{t)[H{t)-i^]Uk{t) with 



Uk{t) 



^~ifi{t)ia^ +a)J^ ^-i^ {,J^+a^ a)t 



(3) 
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^) sin ilt. The explicit calculation of 



where ii{t) 

Hk{t) is reported in Ref. [30|. Here it suffices to note 
that it can be written in the form 



Hk{t) = J2 hl^^ cxp {inilt) . 



(4) 



n— — oo 



In analogy with the standard RWA of quantum optics, 
we obtain an approximate Hamiltonian to describe the 
fcth resonance by neglecting all terms in Hk(t) with os- 
cillatory time-dependence: Hkit) ~ h'^K We therefore 
approximate our time-dependent problem with a series 
of time-independent models, one for each fc > 0. The 
eigenenergies of these time- independent Hamiltonians are 
actually quasienergies [27j . 

We first consider the fundamental fc = case with 
approximate Hamiltonian 

+ 4=(at + a) (,/+ + J_) + 2u: ' :|,(5) 



N 



where S^'^^ = uj—kn/2 and Jq*"^ = ujQ—kQ,/2 are detunings 
in the fcth rotating frame, and J^oix) is the zeroth-order 
Bessel function, here, as in Ref. with operator ar- 
gument. The validity of the approximation Hoit) ~ ft-g"' 
requires 5, w, gjn <^ fl. Omaxd. AQ/17)[30j. In the limit 
— 0, this Hamiltonian recovers the original undriven 
Dicke model. Expanding up to second order in AgJVL, 
we obtain a model similar to the undriven DM 1301 . 



importantly, it opens up new zones of instability in the 
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that can be analyzed in the same way [3|. This anal- 
ysis reveals a second-order QPT along the critical line 
g ~ uj/2 + aj(A(7/r2)^ (a; — ujq), which gives the first 
contribution to the shift of the DM phase boundary ob- 
served in Fig[TJ For the higher resonances, fc > 0, exact 
expressions for /iq'^^ are difficult to obtain. However, it is 
possible to write down explicit expansions up to any finite 
order in Ag/ft. To second order, the k = 1 Hamiltonian 
reads 

The region of validity of this Hamiltonian is w « wo and 
5(^\ 5, [Ag/VLf < oj.uJo. The first two lines of 

this Hamiltonian represent a model similar to the origi- 
nal Dicke model in which energy conserving and noncon- 
serving parts of the interaction have independent cou- 
pling parameters {g and Ag/2, respectively) similar to 
[3l| . At second order in Ag/fi, new effective interactions 
arise: the terms on the third line of Eq. ^ can be in- 
terpreted as an effective dispersive atom-field interaction 
and a dipole-dipole interaction between the atoms. An- 
alytically, this Hamiltonian may again be treated in the 
same way as the undriven DM, from which we observe 
a second-order superradiant transition occurring at the 
critical lines g = -Ag/2 + |(5(i) + u{Ag/nf\. A simi- 
lar analysis for the k = 2 case yields a further second- 
order QPT along the lines at g = 5'^'^^ + ^{Ag/nf and 
g ~ — — ^{Ag/^Y . This analysis can be repeated 
for the all values of fc. 

For small driving, thus. Fig. [T] shows that the original 
DM QPT is joined by a set of new nonequilibrium QPT 
sidebands, the visible number of which (i.e. have signif- 
icant width) increases with increasing driving strength. 
Each of these nonequilibrium QPTs is similar to the orig- 
inal transition in many respects; the transitions are of the 
second-order, mean-field type, with the same critical ex- 
ponents as in undriven DM [3|. An important difference, 
however, is that, whereas the original DM QPT occurs 
for a static coupling g ^ a;/2, the sideband QPTs oc- 
cur for a coupling g ^ <^ w,a;o with driving strength 
also of the same magnitude. In contrast to the static case, 
since the detuning J^*^' can be made arbitrarily small, the 

term Hng = Q' ^^L*!"' ("^ + arising from the square of 
the vector potential can be neglected even for a > 1. 
Therefore, the sideband QPTs are not prohibited by the 
no-go theorem 0,[5[ and should thus be observable in e.g. 
cavity- and circuit- QED experiments [s^ . 

First-order nonequilibrium QPT:- We now revisit the 
full (nonpcrturbative in Ag/^l) fc = resonance Hamil- 
tonian, Eq. ^ . We can investigate the critical properties 
of this model by using the Holstcin-Primakoff representa- 
tion again and by introducing macroscopic displacements 
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FIG. 2: Phase diagram as a function of static coupling 
(g) and driving amplitude (Ag) of effective Hamiltonian 
which describes the driven DM near the fc = resonance. The 
labels indicate the number of local minima of the ground-state 
energy landscape Eq in the corresponding zone. The normal 
phase is the region where there is just a single minimum ( 
X = Y — 0). The superradiant phase is a region with two 
global minima (nonzero order parameter). The boundary be- 
tween these two regions marks a second-order QPT (dotted 
red curve). Outside these regions, the energy surface has an 
odd number > 3 of total minima, sometimes with a single 
global minima, sometimes with two. The boundary between 
these possibilties is a first-order QPT (solid blue line). The 
parameters are cj/Sl = cjo/i^ = 0.05. 

\fxj2 X ia\A^ N 12 Y of the field and atomic modes, re- 
spectively. In the thermodynamic limit, — s- cx), we ob- 
tain a quadratic effective Hamiltonian with leading term 
\XEg(X,Y) with 



Eg{X,Y) = Lu Y^{2 - _ ^XY^/2 - Y'^ 

+.A-+.o(y-lW.(i^). (7) 

The global minima of Eq[X,Y') give the ground-state 

energy of Yl'q^ and either A" or y may be taken as order- 
parameter. In the undriven case, the energy surface 
Eq (X, Y) exhibits a bifurcation from single to double 
minima at the QPT. In the driven case, while retaining 
this bifurcation, Eq (X, y ) also develops multiple addi- 
tional minima, with the general trend that the number 
of minima increases with driving strength (see Fig. [2]). 

Figure [3] shows a sequence of cuts through the en- 
ergy surface for increasing driving strength with static 
coupling g > gc- As is clear from this sequence, at a 
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FIG. 3: Sections of the ground-state energy landscape 
Eg{X,Y) for a fixed value of the static coupling {2g/Q, — 
0.195) and increasing driving amplitudes (Ag) represented 
by dashed straight lines in Fig. [21 This sections correspond 
to the function Eg{X[Y],Y), where X[Y] is a line in the 
order-parameter space crossing all the critical points of the 
ground-state energy. The four panels show (a) the superra- 
diant phase with two global minima, (b) the emergence of a 
local minimum at the origin, (c) multiple minima, but still 
two global minima, and (d) a single global minima at the ori- 
gin plus many further local minima. Panels (c) and (d) are 
separated by the first-order phase transition. The parameters 
are uj/il = ujo/fl = 0.05. 

certain value of Ag, the potential switches from having 
two global minima to just a single one at the origin. At 
the boundary between these behaviors, the potential ex- 
hibits three global minima where normal-phase-like (with 
X = 0) and superradiant-like (with |X|, |F| > 0) so- 
lutions coexist. The system is therefore seen to exhibit 
a first-order QPT [3^ as a function of driving strength, 



the location of which is shown in Fig. [2l Furthermore, 
we interpret the local minima of E{X, Y) as metastable 
phases of h^^K These metastable states are related to 
the phenomenon of parametric stabilization (23 |. Since 
these are separated from the global minima by macro- 
scopic displacements, we expect transitions to the actual 
ground-state to be suppressed, such that the correspond- 
ing values of the order-parameters are observable. This 
possibility is reinforced when one recalls that Eq. ([7]) ac- 
tually describes the lowest quasienergy, which does not 
have the same thermodynamic significance as the lowest 
actual energy. 

In summary, then, we have discussed a driven Dicke 
model through the use of a series of effective Hamiltoni- 
ans obtained under a generalized RWA. For weak driving 
the system exhibits a set of QPT side-bands for which 
the no-go theorem is circumvented. At strong driving the 
long-time dynamics of the system in the k = resonance 
case is governed by the effective Hamiltonian /iq"'' , which 
exhibits rich structure with a first-order quantum phase 
transition and metastable states. We anticipate that the 
higher resonances show similar behavior. Our methodol- 
ogy should be generalisable to investigate similar regimes 
for other phase transitions. 

A possible experimental realization of our model could 
be carried out by means of a quantum well waveguide 
structure, which allows one to switch on the light-matter 
interaction within less than one cycle of light An- 
other possible implementation consist in a Bose-Einstein 
condensate coupled to an optical cavity where the effec- 
tive atom-field coupling could be controlled externally by 
varying the intensity of the pump laser as a function of 
time about a static value. 
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Supplementary IMaterial for "Non-equilibrium Quantum Phase Transitions in the 

Dicke Model" 

A. The Hamiltonian in the rotating frame 

We study the general single-mode driven Dicke model: 



Hit) = coa'' a + too Jz + "^{a^ + a) {J+ + J-), 

V N 



(8) 



and the case of a monochromatic modulation of the atom-field coupling g{t) = g + Ag cos fit. We perform a canonical 
transformation of the Hamiltonian Eq. ([5]) into a rotating frame using 



Ukit) = exp 



2 / AgsinOt\ 



cxp 



(9) 



In this frame, we study the quantum dynamics using the Hamiltonian Hk{t) = ul{t){H{t) - i§^)Uk{t). The explicit 
form of Hk{t) is 



N 

Lo ( IS.g sin fii 



2 / A(7siur2t\ / i ikvit iksit 

■= a^e— + ae — 

N \ ^ /V 

kilt \ ( ikilt -iknt 

ae 2 ) ( J, e 2 + j g 2 



kn/2j Jz 



a'e 



iknt \ ( ikilt -iknt 

2 j ( J+e 2 -I- J_e 2 



( AqsinOtX / iknt -iknt 



.^0 . 
-t — sm 
2 



2 / Agsinr2t\ / + iknt iknt 
^(^^](a^e^+ae-^ 



where 5^^^ ^ uj — k^l/2. We then rewrite this equation using the operator identities 

oo 

cos{d{t)smnt) = Jo{d{t)) + 2^J2n{d{t))cos{2nnt), 

n=l 

oo 

sm{d{t)smnt) = 2^J2n+i{d{t))sm{{2n + l)nt) , 



iki'2t -iki'lt 

Ji e 2 - J_e 2 



(10) 



(11) 



n=0 
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where J7m(- ) denotes the Bessel function of integer order m and the time-dependent operator argument is given by 
0{t) = ('^) ("^^6^^ + ae~^^y Considering the power series expansion of the Bessel functions, it is possible 
to write the Hamiltonian Eq. (fTO|) in the general form 



oo 

Hk{t)^ hl^^ exp{innt). (12) 

n— — oo 

In the rotating frame we are able to describe explicitly multiphoton processes that arise as a consequence of the 
interplay between the quantum cavity field and the classical external driving [2[ . 

To perform the generalized RWA Q-Q , we neglect the contributions due to virtual processes that ocurr faster than 

the characteristic time scales of the system. Therefore, we consider only the the zero frequency component h'^^^ of 
the Hamiltonian Eq. For example, in the case of the k = resonance, a general condition of RWA validity 

in the strong driving regime should be satisfied for all \n\ > 0: g, ujo\\Jn{0{0))\\, wj] J„((!)(0))|| <C f^, Ag, where ||-|| 
is an operator norm. In contrast, in the weak driving regime A.g/D, <C 1, we require g, lu, loq <ti The norm of 
the Bessel function is always bounded, therefore we can summarize the conditions for weak and strong driving as 
g, w, ujq <^ n^fl ■ max{l, Ag/n). 

The full fundamental fc = effective Hamiltonian drived with this approach reads: 



<5Wata + C^of;^(at + a)') J. 



--^(at+a)(J, + J_) + 2.(^)'|. (13) 



In the weak driving regime, the approximation to this Hamiltonian li'^^ up to second order in Agj^ is 



N 



To second order, the fc = 1 effective Hamiltonian reads 



fAgV ^ J (AgVj^J+ 



and similarly for the fc = 2 case 



n J N \n J N 



I. THE THERMODYNAMIC LIMIT 

To address the expansion of the effective Hamiltonian h^^^ we begin by considering the Holstein-Primakoff repre- 



sentation of the angular momentum algebra @ 



= b^b - N/2, J- = Vn - 5t6 6, J+ = bWN -b^b. 



(17) 
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In this representation, the driven Dicke Hamiltonian is written in terms of two bosonic modes a and h. We introduce 



macroscopic displacements X^jN Y^J N/2 of order \/N for the bosonic a and for the atomic b modes, respectively. 
These macroscopic displacements are defined by 

a = c±X^-, b = dTY^-. (18) 

In order to describe the QPT, we derive effective quadratic bosonic Hamiltonians for both the normal phase and the 
superradiant phase 0- In the thermodynamic limit N — > oo, we perform a series expansion of the Hamiltonian h^^'' 
m powers of 

hl,'^ « /ig) (c, d,c\d^) + ^ hf (c, d,c\d^) + ^ Eg{X,Y), (19) 

where K^q' is the desired effective quadratic bosonic Hamiltonian (depending on the choice of the macroscopic dis- 
placements Eq. (lisp ), hf contains linear bosonic terms, and Eij{X^Y) is the lowest quasienergy. 
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